Engineering a robust quantum spin Hall state in graphene via adatom deposition 
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The 2007 discovery of quantized conductance in HgTe quantum wells delivered the field of topological insu- 
lators (TIs) its first experimental confirmation. While many three-dimensional TIs have since been identified, 
HgTe remains the only known two-dimensional system in this class. Difficulty fabricating HgTe quantum wells 
has, moreover, hampered their widespread use. With the goal of breaking this logjam we provide a blueprint for 
stabilizing a robust TI state in a more readily available two-dimensional material — graphene. Using symmetry 
arguments, density functional theory, and tight-binding simulations, we predict that graphene endowed with 
certain heavy adatoms realizes a TI with substantial band gap. For indium and thallium, our most promising 
adatom candidates, a modest 6% coverage produces an estimated gap near 80K and 240K, respectively, which 
should be detectable in transport or spectroscopic measurements. Engineering such a robust topological phase 
in graphene could pave the way for a new generation of devices for spintronics, ultra-low-dissipation electronics 
and quantum information processing. 

PACS numbers: 



I. INTRODUCTION 

When certain physical properties of a system depend on 
global topology — and not on local details, such as disorder — 
then the system is said to realize a topological phase. By 
virtue of their universality, rooted in this independence on 
microscopic details, topological phases are both paradigmat- 
ically interesting and thought to possess seeds of impor- 
tant future technologies. Until rather recently, experimen- 
tal studies of topological phases were confined to the quan- 
tum Hall effect, seen in high-quality two-dimensional elec- 
tron systems subjected to strong magnetic fields. The ad- 
vent of topological insulators (TIs), a class of two- and three- 
dimensional non-magnetic crystalline solids enjoying strong 
spin-orbit couplin^^, opened a fascinating new chapter in 
the field. 

Two-dimensional TIs are commonly referred to as quan- 
tum spin Hall (QSH) states' ^ due to their similarity with in- 
teger quantum Hall liquids. A defining signature of a QSH 
state is the existence of gapless edge states which are pro- 
tected from elastic backscatteiing and localization by time- 
reversal symmetry (T). This feature, combined with the spin- 
filtered nature of the edge modes, renders such systems tech- 
nologically very promising. Furthermore, when in contact 
with a superconductor the edges are predicted to host Ma- 
jorana fermions**, which play an important role in topologi- 
cal quantum computation schemes^. The QSH state was first 
predicted to arise in graphene*", though later it was realized 
that, due to carbon's small atomic number, spin-orbit cou- 
pling is too weak to produce an observable effect under re- 
alistic conditions. (First-principles calculati on J^El on pristine 
graphene all predict sub-Kelvin gaps for this phase). The QSH 
state was subsequently predicted in several ot her m aterials, in- 
cluding HgTe and InAs/GaSb quantum well^^^EH and bilayer 
bismuth'^. While landmark experiments indeed observed this 
phase in HgTe'^, experimental activity on the QSH effect has 
remained limited by notorious practical difficulties associated 
with this material. 



Given the comparative ease with which graphene can be 
fabricated, it would be highly desirable to artificially enhance 
graphene's spin-orbit coupling strength to elevate the gap pro- 
tecting the QSH state to observable levels. A practical means 
of doing so would pave the way to a new generation of QSH 
experiments accessible to a broad spectrum of researchers 
worldwide. Previous work has explored the possibility that 
strong interactionJi^GS or radiation'^ might drive graphene 
into a robust topological phase. Our goal in this manuscript 
is to explore a practical new route to this end. Specifically, we 
will ask whether graphene can 'inherit' strong spin-orbit cou- 
pling from a dilute concentration of heavy adatoms (whose 
innate spin-orbit coupling strength can be on the electron-Volt 
scale) deposited randomly into the honeycomb lattice. 

The basic principle underlying our proposal can be under- 
stood by considering processes in which an electron from 
graphene tunnels onto an adatom — whereupon it 'feels' enor- 
mous spin-orbit coupling — and then returns to the graphene 
sheet. Such processes in effect locally enhance graphene's 
spin-orbit strength manyfold. A viable proposal must, how- 
ever, address numerous competing effects. For instance, 
adatoms often form local magnetic moments-" which poten- 
tially spoil the time-reversal symmetry protecting the QSH 
effect. Moreover, adatoms generically mediate both intrin- 
sic and Rashba spin-orbit coupling. The latter is believed to 
be detrimental to the QSH phase^, and previous work has in- 
deed established that certain kinds of adatoms do generate 
substantial Rashba coupling in graphene that typically over- 
whelms the intrinsic contributiorf^"*^^'. (As Ref. '2T showed, 
however, magnetic adatoms inducing strong Rashba coupling 
can induce an interesting 'quantum anomalous Hall' state in 
graphene.) The adatoms may also favor competing, ordinary 
insulating states depending on their precise locations in the 
lattice. And finally, since spin-orbit coupling is generated 
non-uniformly in graphene, the stabilization of a QSH phase 
even in an otherwise ideal situation is unclear a priori. 

After an extensive search using tight-binding and first- 
principles analyses, we have found that two elements, indium 
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FIG. 1: Adatoms in graphene. (a) Depending on the element, 
adatoms favor either the high- symmetry 'bridge' (B), 'hollow' (H), 
or 'top' (T) position in the graphene sheet, (b) Detailed view of an 
H position adatom, which is best suited for inducing the intrinsic 
spin-orbit coupling necessary for stabilizing the topological phase. 
The desired spin-orbit terms mediated by the adatom are illustrated 
in (c). Red and yellow bonds represent the induced second-neighbor 
imaginary hopping, whose sign is indicated by the arrows for spin up 
electrons. For spin down electrons the arrows are reversed. 



and thallium, are capable of stabilizing a robust QSH state in 
graphene. Neither element forms a magnetic moment, and 
although they do generate significant Rashba coupling, for 
symmetry reasons this remarkably does not suppress the QSH 
state. We find that gaps many orders of magnitude larger than 
that predicted in pure graphene can form even with coverages 
of only a few percent; for example, at 6% coverage indium 
yields a gap on the order of lOOK, while for thallium the 
gap approaches room temperature. These predictions revive 
graphene as a viable QSH candidate, and can be verified by 
probing the gap and associated edge states using spectroscopic 
and conductance measurements. 



II. PHYSICS OF A SINGLE HEAVY ADATOM 



spin indices, Ht and Hso explicitly read 

Ht = ^ (4cr' + /i.e.) (1) 

(rr') 

Hso = Aso ^ {iVrr'cls^Cr' + h.c). (2) 
((rr'» 

Here and below, spin indices are implicitly summed when- 
ever suppressed. In Eq. (j2]i the sum runs over second-nearest- 
neighbor lattice sites, and v^^' are signs that equal +1 if an 
electron hops in the direction of the arrows in Fig. [TJc) and 
— 1 otherwise. Thus H^o describes 'chiral' spin-dependent 
second-neighbor electron hopping. When Ago = 0, the band 
structure exhibits the familiar gapless Dirac cones centered on 
momenta ±Q, resulting in semimetallic behavior Turning on 
Ago 7^ generates an energy gap* A ~ 6\/3|Aso| at the Dirac 
points, transforming the system into a (very fragile^f''') QSH 
insulator Importantly, if miiTor symmetry with respect to the 
graphene plane is broken, then Rashba coupling — which in- 
volves spin flips and thus breaks the U(l) spin symmetry en- 
joyed by H — will also be present*. Rashba coupling com- 
petes with the intrinsic spin-orbit term in pure graphene, and 
beyond a critical value closes the gap and destroys the QSH 
state. 

If heavy adatoms are to stabilize a more robust QSH phase 
in graphene, then at a minimum they should be non-magnetic 
(to preserve 7") and modify the physics near the Dirac points 
primarily by inducing large intrinsic spin-orbit coupling. The 
latter criterion leads us to focus on elements favoring the 'hol- 
low' (H) position in the graphene sheet indicated in Fig. [TJa). 
Compared to the 'bridge' (B) and 'top' (T) positions, adatoms 
in the H position can most effectively mediate the spin- 
dependent second-neighbor hoppings present in the Kane- 
Mele model, while simultaneously avoiding larger competing 
effects such as local sublattice symmetry breaking generated 
in the T case. 

Since H-position adatoms generically reside on one side of 
the graphene sheet, they will clearly mediate Rashba spin- 
orbit coupling as well, leading to a potentially delicate com- 
petition. If the adatoms' outer-shell electrons derive from p- 
orbitals, however, the induced intrinsic spin-orbit terms al- 
ways dominate over the induced Rashba interactions. One can 
establish this key result by studying graphene with a single 
adatom of this type. To model this setup, we employ operators 
dma for the adatom states, with m = 0, ±1 and a =t, i label- 
ing the orbital and spin angular momentum quantum numbers. 
The coupling of these orbitals to graphene is conveniently ex- 
pressed in terms of the following operators: 



To set the stage, let us first briefly review the Kane-Mele 
model* describing pure, undoped graphene with spin-orbit 
coupling. The Hamiltonian can be expressed as iJKM = 
Ht + Hso, where Ht describes the usual nearest-neighbor hop- 
ping and Hso encodes intrinsic spin-orbit coupling. In terms 
of operators cj^ that add electrons with spin a to site r of 
the honeycomb lattice and Pauli matrices s^-V''^ that act on the 



where the sum runs over the six sites surrounding the adatom 
shown in Fig. [ijb). Guided by symmetry (see Appendix 
A), we consider the following minimal Hamiltonian for this 
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single-adatom problem: 

H = Hg + Ha + H, (4) 

6 

Hg = Ht - Sfi^clxrj (5) 

m=0,±l 

+ V2A'^^{dls-d^i+dls+di + h.c.) (6) 
He = - J2 {t\m\Cid„, + h.c.), (7) 

m=0,±l 

with — (s^ ± is^)/2. Here represents the nearest- 
neighbor hopping model of graphene supplemented by a 
chemical potential 6fi for the six sites surrounding the adatom. 
Physically, a non-zero dj.! leads to an excess electron density at 
those sites, screening any net charge from the adatom. Crystal 
field effects and spin-orbit coupling split the adatom p-orbitals 
through the e|„j| and AscAg^ terms in Ha- Finally, He al- 
lows electrons to tunnel between the adatom and its neigh- 
boring carbon sites. All couplings in H are real except <i, 
which is pure imaginary. Note also that we have ignored the 
exceedingly weak spin-orbit terms that couple the electrons 
in graphene directly, as well as Coulomb interactions for the 
adatom. Exclusion of the latter terms is justified for the non- 
magnetic adatoms we will consider below. 

Recalling that intrinsic spin-orbit coupling in graphene con- 
serves the 5*^ spin component while Rashba coupling does 
not, one can readily understand how the adatom mediates 
both kinds of interactions by viewing the tunnelings in He 
perturbatively. For example, when an electron hops onto the 
m — orbital via t^, flips its spin through the A^^ coupling, 
and then hops back to the honeycomb lattice via ti, spin-orbit 
coupling of Rashba type is locally generated in graphene. In 
sharp contrast to the situation for the Kane-Mele model, how- 
ever, Rashba terms mediated in this fashion are irrelevant for 
the low-energy physics. This crucial feature can be under- 
stood by Fourier transforming Cqq; remarkably, the compo- 
nents with Dirac -point momenta ±Q vanish identically. More 
physically, electrons near the Dirac points interfere destruc- 
tively when hopping onto the m = adatom orbital. (See 
Appendix A and D for a complementary perspectives.) We 
stress that this argument holds only for p-orbitals. If the rel- 
evant adatom orbitals carry higher angular momentum, spin- 
flip processes which do affect the Dirac points will generi- 
cally appear. The induced Rashba terms may still be subdom- 
inant, though whether this is the case depends on details of the 
Hamiltonian unlike the situation for p-orbitals. 

In contrast, 5^ -conserving events whereby an electron hops 
both onto and then off of an m = ±1 orbital via ti locally 
mediate the desired intrinsic spin-orbit interactions. No ob- 
struction exists for tunneling onto the m = ±1 orbitals via 
ti, so these processes can effectively modify the physics near 
the Dirac points. It is important to note, however, that ti 
hopping mediates additional couplings as well, so the dom- 
inance of these spin-orbit terms remains unclear. Neverthe- 
less, given the high symmetry preserved by the H-site adatom. 



none of the conventional broken-symmetry gapped phases of 
graphene — such as charge-density wave or 'Kekule' orders — 
are obviously favored here, so it is reasonable to expect these 
additional terms to play a relatively minor role. We explicitly 
confirm this intuition in the multi-adatom situation, to which 
we now turn. 



III. PERIODIC ADATOM CONFIGURATIONS 

As a first step in understanding the multi-adatom case, we 
examine a periodic system with one adatom residing in a large 
N X N supercell. This situation allows us to utilize den- 
sity functional theory (DFT) to ascertain suitable heavy ele- 
ments and obtain a quantitative understanding of their effects 
on graphene. (See Appendix B for computational details.) To 
ensure large spin-orbit coupling, we focused on elements in 
rows five and six of the periodic table, including In, Sn, Sb, 
Te, I, La, Hf, Pt, Au, Hg, TI, Pb and Bi. For each element, we 
calculated the total energy in the three adsorption geometries 
shown in Fig. 1(a) along with the adatom's spin moment. Our 
calculations reveal that two elements — indium (atomic num- 
ber Z — 49) and thallium (Z = 81) — satisfy our criteria of 
both favoring the high-symmetry H position and being non- 
magnetic. Furthermore, both elements exhibit partially filled 
p-shells, ensuring that the Rashba coupling they mediate in 
graphene is benign at the Dirac points. 

It is instructive to first examine the electronic properties of 
indium on graphene in a 4 x 4 supercell, without spin-orbit 
coupling. As Fig. [2|a) illustrates, the Dirac cones character- 
istic of pure graphene indeed remain massless — despite the 
reduced translation symmetry, conventional gapped phases 
are not stabilized here, consistent with the intuition devel- 
oped in the single-adatom case above. Indium does, however, 
electron-dope graphene and shifts the Fermi level Ep 0.95 
eV above the Dirac points. From the adatom's local density 
of states (LDOS) displayed in Fig. |2|a), one can see that in- 
dium's bp orbitals lie almost entirely above Ep, implying that 
the 5p electron in neutral indium nearly completely transfers 
to graphene. (The charge of an indium adatom is +0.8e from 
the Bader charge division scheme.) Note that the relatively 
diffuse pz LDOS indicates that this orbital hybridizes more 
strongly with graphene compared to the p^^y orbitals. Re- 
placing indium with thallium, again without spin-orbit cou- 
pling, leads to the band structure and LDOS shown in Fig. 
|2|d). Clearly the electronic structure is modified very little by 
this substitution; importantly, the Dirac cones remain mass- 
less with thallium as well. 

Thus any gap opening at the Dirac points must originate 
from spin-orbit coupling. Fig.|2|b) displays the band structure 
and LDOS for spin-orbit-coupled indium on graphene. Note 
the sizable spin-orbit splitting in the LDOS for the p^.^j or- 
bitals. More remarkably, a gap Ago ~ 7meV now appears 
at the Dirac points, which already exceeds the spin-orbit- 
induced gap in pure graphen^^^ by several orders of mag- 
nitude. The analogous results for thallium — whose atomic 
mass is nearly twice that of indium — are still more striking. 
As Fig. |2|e) illustrates, p-orbital splittings of order leV are 
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FIG. 2: Band structure and the adatom local density of states (LDOS). All data correspond to one adatom in a 4 x 4 supercell, with 
the upper row corresponding to indium and the lower row corresponding to thallium. The left panels in (a) and (d) correspond to the band 
structure and LDOS computed using DFT without spin-orbit coupling. The horizontal dashed lines indicate the Fermi level {Ef), which 
shifts due to electron-doping from the adatoms. Insets zoom in on the band structure near the K point within an energy range —35 to 35meV, 
showing that without spin-orbit interactions neither indium nor thallium open a gap at the Dirac points. The central panels in (b) and (e) are 
the corresponding DFT results including spin-orbit coupling. Remarkably, in the indium case a gap of 7meV opens at the Dirac point, while 
with thallium the gap is larger still at 21meV. Finally, the right panels in (c) and (f) were obtained using the tight-binding model described in 
the text. 



now evident in the LDOS, and a gap Ago ~ 21meV opens at 
the Dirac points. We emphasize that these results apply for 
adatom coverages of only 6.25%. To explore the dependence 
of Ago on the adatom coverage, we additionally investigated 
systems with one adatom in 5 x 5, 7 x 7 and 10 x 10 su- 
percells. (For N x N supercells with N a multiple of 3, the 
Dirac points reside at zero momentum and can thus hybridize 
and gap out even without spin-orbit coupling. We thus ignore 
such geometries.) The values of Ago along with the Fermi 
level Ep computed for the coverages we studied appear in 
Fig. [3] Quite naturally, the gap decreases as one reduces the 
coverage, as does the Fermi level since fewer electrons are do- 
nated to graphene at lower adatom concentrations. It is worth 
highlighting that the gap decreases sublinearly and that a siz- 
able Ago ~ lOmeV remains even with a mere 2.04% thalhum 
concentration. 

Since spin-orbit coupling clearly underlies the gap onset, 
it is tempting to associate the insulating regime opened by 
the adatoms with a QSH state. To verify this expectation, 
we appeal to tight-binding simulations of the 4x4 super- 
cell geometry. We model the system by a Hamiltonian i?4x4 
consisting of uniform hopping for graphene (including weak 
second- and third-neighbor tunneling) supplemented by the 



additional local terms in Eq. Q at each adatom. Figure [2|c) 
displays the tight-binding band structure and adatom LDOS 
for indium with the following parameters: first, second, and 
third-neighbor graphene hoppings t ~ 2.82eV, t' = 0.22eV, 
and t" = 0.2eV; 6^ = 0.5eV; eo = 2.5eV, ei = 1.8eV; 
to = 2eV, iti = 0.95eV; and A^o = = O.leV. The 
corresponding data for thallium appear in Fig. |2]^f); param- 
eters are the same as for indium, except for ei ~ 1.9eV 
and Ago = A^^ = 0.31eV. The good agreement with DFT 
provides strong evidence that the underlying physics is well- 
captured by our tight-binding model. With this in hand, one 
can then demonstrate the topological nature of the insulating 
regime by showing that the Dirac point gap remains non-zero 
upon adiabatically deforming the Hamiltonian to the Kane- 
Mele model. This can be achieved by defining a new Hamil- 
tonian H{X) = (1 — A)iJ4x4 + Ai?KM that interpolates be- 
tween our adatom model at A = and the Kane-Mele model 
at A = 1. For either indium or thallium parameters, the Dirac 
point gap never closes as A varies from to 1 (see Appendix C 
for details). Thus either element can indeed stabilize a robust 
QSH state in graphene, at least when periodically aiTanged. 
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FIG. 3: Coverage Effects. Spin-orbit-induced gap Ago opened at 
the Dirac points and Fermi level Ef (measured relative to the cen- 
ter of the gap) for different indium and thallium adatom coverages. 
The open and filled symbols represent data for indium and thallium, 
respectively. 



IV. RANDOMLY DISTRIBUTED ADATOMS 

In an experiment adatoms will likely occupy random Hi- 
site locations in a graphene sheet that is of course not pris- 
tine. It is therefore important to understand the impact of such 
randomness and disorder on the stability of the topological 
phase that was found above in pure graphene with periodic 
adatoms. To address general features of the problem in a way 
that reduces its computational complexity, we will describe 
the disordered, random-adatom system using the following 
graphene-only Hamiltonian, 



(8) 



Here Sfir represents a random on-site potential arising, e.g., 
from the substrate (and not the adatoms). In the second term, 
/ labels the random plaquettes occupied by adatoms and 

dHj — — (5/1^^ c|!Cr + Aso {iVrr'cls^Cr' + h.C.) 

re/ ((rr'))e/ 

+iXR ^ 4(S X drrO^Cr' (9) 
T.r'el 

with drr' = r — r'. The first term in 5Hi describes the chem- 
ical potential that screens charge from the adatoms, while the 
last two terms capture the local intrinsic (Xso) and Rashba 
(A r) spin-orbit couplings induced by electrons hopping from 
graphene to an adatom and back, as discussed in Sec. 11. For- 
mally, these spin-orbit terms can be derived by integrating out 
perturbatively the adatom degrees of freedom in Hamiltonian 
(5]|7i; see Appendix D. It is important to notice that unlike 
the conventional nearest-neighbor Rashba term considered in 
Ref. 6 for pristine graphene, the adatom-generated Rashba 
term connects all sites in the hexagon. Such a 'hexagon 
Rashba term' has the property of vanishing at the Dirac 
points as akeady argued in Sec. 11. The adatom also induces 
other symmetry-allowed terms, such as further-neighbor spin- 
independent hoppings, which we disregard because they are 



either weak or do not lead to qualitative changes in the results 
reported below. 

To study the robustness of the QSH phase in the presence of 
random adatoms and disorder we probe for the signature gap- 
less edge states by simulating the classic two-terminal con- 
ductance measurement in the geometry of Fig. |4ja). Within 
the Landauer-Biittiker formalisrrl^'^we employ the recursive 
Green's function method^'' to evaluate the conductance G of 
a length-i, width- graphene strip as we vary the Fermi en- 
ergy. Consider first the simplest random-adatom setup with 
(5/ir — (corresponding to an otherwise clean graphene 
sheet), Sfi = 0, and A^ = 0. Figure |4jb) illustrates the con- 
ductance for this case at several adatom concentrations n^, and 
model parameters indicated in the caption. A 2e^/h plateau 
with width proportional to n,; clearly emerges, strongly sug- 
gesting the onset of a bulk mobility gap and quantized ballistic 
conduction by edge states. This picture is corroborated by Fig. 
|4|c), which displays the current distribution for Ep ~ Q.15eV 
across a smaller system size (chosen for clarity). 

These results establish that in principle adatoms need not 
be periodic to stabilize a QSH phase, though the parameters 
employed were unrealistic. To make contact with our DFT 
results, we also considered a system with Aso = 0.02< and 
(5/i = O.lt, which would yield a similar bulk mobility gap to 
that seen in our thallium simulations if the adatom coverages 
were the same. While finite-size effects prevent us from study- 
ing the low-adatom coverages considered in Sec. Ill, Fig. 4(d) 
shows that for Ui ~ 0.15 a robust conductance plateau in- 
deed persists for these more realistic parameter values. (At 
much lower coverages the width of the edge states approaches 
the system sizes we are able to simulate.) We have also stud- 
ied the effects of the hexagon Rashba term and residual on- 
site potential disorder on the stability of the QSH phase (see 
Appendix E for details). Even a relatively large Xr ^ 2Aso 
has a weak effect on the width of the conductance plateau, 
as expected given our analysis in Sec. 11. With uncorrected 
on-site disorder (5/ir, the topological phase is also remarkably 
robust — (5/ir can fluctuate on the scale of t while degrading 
the plateau only marginally. The topological phase is more 
sensitive, however, to correlated disorder which is likely more 
relevant experimentally. In this case the conductance plateau 
survives only when Sfi^ varies on a scale smaller roughly than 
the mobility gap for the clean case. To observe the QSH phase 
experimentally, the chemical potential should therefore fluc- 
tuate on scales smaller than the ^ lOmeV gaps we predicted 
here. While this might be difficult to achieve with graphene 
on standard Si02 substrate, recent experiments using hexago- 
nal boron nitride substrates^^ show disorder energy scales as 
low as 15K, which should be sufficient to observe the QSH 
gap- 



V. RELATION TO PREVIOUS WORK 

Since numerous recent studies have explored the possibility 
of employing adatoms to control graphene's electronic proper- 
ties, here we briefly comment on the distinctions between our 
proposal and related earlier works. Castro Neto and Guine^^ 
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FIG. 4: Transport, (a) Configuration used to determine the two-terminal conductance. Semi-infinite clean graphene leads are connected to 
the sample, which has adatoms distributed randomly across it. The length L of the system is set by the number of columns (indicated by dotted 
lines). The width W refers to the number of sites along either edge of the figure, corresponding to the number of unit cells of an armchair 
nanoribbon required to create an individual column, (b) Conductance G as a function of the Fermi energy Ef, averaged over 40 independent 
random adatom realizations at different concentrations = 0.1,0.2,0.3 for a system of size = 80 and L = 40 with Ago ~ O.lt. (c) 
Current distribution across a sample of size = 40 and L = 20 at rii = 0.2, Ago = 0-U and Ef ~ 0.15eV. The magnitude of the current 
is represented by both the arrow size and color (d) Conductance for the largest simulated system size using realistic parameters for thallium 
adatoms (Aso ~ 0.02t and Sjj. — O.lt) estimated from DFT data. Here W — 200, L = 100, and the coverage is Ui = 0.15. In all panels we 
take t = 2.7eV. 



previously argued that graphene's spin-orbit strength can be 
greatly enhanced by adatom deposition, though the physical 
mechanism and nature of induced spin-orbit coupling is quite 
different from what we considered here. The enhancement 
discussed in Ref. 21 stemmed from structural changes in the 
honeycomb lattice induced by a T-position adatom (such as 
hydrogen whose innate spin-orbit coupling is weak). We ex- 
plicitly neglected such contributions, as our DFT simulations 
demonstrated that the adatoms we considered only weakly af- 
fect the carbon-carbon bonds. Perhaps even more importantly, 
the effect was discussed solely in the context of enhanced 
Rashba spin-orbit interactions, which do not stabilize a QSH 
state in graphene. 

More recently, Abdelouahed ef a/.'^examined the influence 
of heavy Au adatoms on graphene's electronic structure using 
first-principles calculations. These authors considered dense 
coverages, with one T-position adatom residing above each 
carbon site on graphene's A sublattice, and found substantial 
enhancement of graphene's Rashba coupling to values of or- 
der lOmeV. Again, this type of induced spin-orbit coupling is 
not suitable for driving a QSH state, and in any case this ge- 
ometry is very likely to produce only a topologically trivial 
band insulator due to explicit breaking of sublattice symme- 



try by the adatoms. (Lighter Ni adatoms were also studied, 
but were found to modify graphene's spin-orbit coupling only 
very weakly.) 

Closest in spirit to our study is the very interesting work of 
Qiao ef al^, who considered the possibility that adatoms can 
stabilize another topological phase in graphene — a quantum 
anomalous Hall state. Unlike the QSH phase, this state breaks 
time-reversal symmetry and exhibits chiral edge modes lead- 
ing to a non-trivial Hall conductivity as in the integer quantum 
Hall effect. Thus the requirements for stabilizing a quantum 
anomalous Hall phase in graphene are quite different. Qiao ef 
al. showed that this can be accomplished by a combination of 
strong Rashba spin-orbit coupling and an exchange-induced 
Zeeman field. They further demonstrated via first-principles 
calculations that 6.25% coverage of periodically arranged iron 
adatoms (which form magnetic moments) provide both ingre- 
dients and generate a quantum anomalous Hall phase pro- 
tected by a 5.5meV gap. Note that while iron favors the H 
site like indium and thallium, its outer-shell electrons derive 
from d rather than p orbitals. Thus our argument for the irrel- 
evance of Rashba spin-orbit coupling at the Dirac points does 
not apply to iron adatoms, consistent with the numerical find- 
ings of Ref. ,23, Although the gap produced by iron is smaller 
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than that produced by indium and especially thallium at the 
same coverage, a promising feature of their proposal is that 
iron adatoms do not alter the Fermi level. 

We also wish to point out a connection of our work to 
'topological Anderson insulators' — systems which are non- 
topological in the clean limit but acquire non-trivial topology 
when disordered. Such states were first predicted in Ref. |281 
and explored further in several subsequent studieJ^lEl ^hg 
robust QSH state that can emerge upon randomly depositing 
dilute adatoms onto pure graphene (which by itself is a non- 
topological semimetal) can be viewed as an example of a topo- 
logical Anderson insulator. Note, however, that this catego- 
rization is somewhat loose, as this phase in fact smoothly con- 
nects to the topological phase of the pure Kane-Mele modeP, 
which has no disorder (see also Ref. i32l . Nevertheless, this 
viewpoint is useful in that it highlights the remarkable fact 
that disorder, with suitable properties, can provide a practical 
avenue to generating topological phases. 



VI. DISCUSSION 

In this paper we have demonstrated that dilute heavy 
adatoms — indium and thallium in particular — are capable of 
stabilizing a robust QSH state in graphene, with a band gap 
exceeding that of pure graphene by many orders of magni- 
tude. Aside from fabrication ease, another virtue graphene 
possesses as a potential QSH system is the large number of ex- 
perimental probes available in this material because the elec- 
trons are not buried in a heterostructure. Angle resolved pho- 
toemission spectroscopy (ARPES), for instance, actually ben- 
efits from electron doping by the adatoms, and thus could 
be employed to detect the onset of the Dirac point gap at 
arbitrary coverages. Scanning tunneling microscopy (STM) 
would be similarly well-suited for probing the bulk gap, the 
spatial structure of the LDOS around an adatom, and even the 
topologically protected gapless edge modes. 

Detecting the edge states via transport will be more chal- 
lenging, since this would require repositioning the Fermi level 
inside of the bulk gap. Nevertheless, the technology for dop- 
ing graphene by the requisite amount does exist (see, e.g., 
Refs. 33-37). Although the most commonly used ion-liquid 
techniques-'^ •'^ work best near room remperature, recently de- 
veloped approaches employing ferroelectric substrates^^ and 
solid polymer electrolytes^^ show comparable results and re- 
main applicable down to cryogenic temperatures. The requi- 
site doping may also be achievable by introducing additional 
high-electronegativity adatoms that absorb electrons. A quan- 
tum spin Hall phase in an appropriately doped graphene sam- 
ple can be detected through the classic two-terminal conduc- 
tance measurement akin to the historic observation performed 
in HgTe quantum wellJ*^. A more complete characterization 
of the topological phase can be obtained through non-local 
edge state-''^ and magnetotransport measurements which pro- 
vide unambiguous evidence for topologically protected heli- 
cal edge states-'^'*". 

Fulfilling Kane and Mele's original vision by artificially 
turning graphene into a strong spin-orbit system could have 



remarkable technological implications. Among these, the 
possibility of employing the protected edge states for real- 
izing topological superconductivity, Majorana fermions, and 
related phenomena such as non-Abelian statistics and 'frac- 
tional Josephson effects' are particularly tantalizing"*. In- 
terestingly, one of the primary requirements for observing 
these phenomena — generating a superconducting proximity 
effect in graphene — has already been achieved by numerous 
groups'*' These exciting possibilities, along with potential 
spintronics and low-dissipation devices, provide strong moti- 
vation for pursuing the line of attack we introduced for stabi- 
lizing a robust QSH state in graphene. 
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Appendix A: Discussion of tlie effective Hamiltonian for 
grapliene witli a single adatom 

In this section we will consider a sheet of graphene located 
in the (x, y) plane, with a single adatom sitting at the H posi- 
tion (recall Fig. 1 of the main text). An important feature of 
such an adatom is the high symmetry that it preserves com- 
pared to the cases where the adatom resides at the B or T po- 
sitions. In particular, the following symmetries remain here: 
7r/3 rotation {Rtt/z) and x-reflection {Rx) about the adatom 
site, as well as time-reversal (T). Note that since the adatom 
is displaced vertically out of the {x, y) plane, mirror symmetry 
with respect to the graphene plane {Rz) is broken. It will nev- 
ertheless be useful below to understand how operators trans- 
form under R^, since this will give us insight as to which cou- 
plings mediate intrinsic versus Rashba spin-orbit coupling. 

We will continue to assume that p-orbitals form the active 
states for the adatom. Also as in the main text, we employ 
operators dj^j^, which add electrons with spin a and orbital an- 
gular momentum quantum number rn = 0, ±1 to the adatom, 
and operators cj^ which add electrons with spin a to site r 
of graphene's honeycomb lattice. These operators transform 
under ii!^/3 according to 

i?^/3 : Cr e*^"'cr' (Al) 

i?V3 ■ ^™ ^ e'^^'e'T^d™ (A2) 

where r' corresponds to r rotated by 7r/3. Note that here and 
below, both Cr and d,n are considered to be two-component 
objects, with the upper and lower components corresponding 
to spin up and down, respectively. Similarly, under R^ we 
have 

Rx ■■ C(^x,y) e*"5""c(_^,y) (A3) 

Rx- dm ^ e'^'^'d-™. (A4) 
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Time-reversal, which is antiunitary, sends 

T: Cr ^ e'^'\r (A5) 

T: d,„ ^ (-l)'"e'i^''d_,„. (A6) 

Note that the (—1)™ factor above arises because one con- 
ventionally defines spherical harmonics such that [Y^™]* = 
(— l)™y^^'". Finally, the operators transform under as 

R^ : Cr ^ -e^^^'cr (A7) 

R. : dm ^ (-l)i-'"e'5^^d„ (A8) 

In the last equation, the factor of (— reflects properties 
of spherical harmonics under z — > — z. 

Our aim now is to construct a minimal tight-binding model 
for graphene with a single H-position adatom. Following the 
main text, we write the Hamiltonian as H — Hg + Ha + He, 
where Hg involves only the carbon degrees of freedom. Ha 
involves only the adatom degrees of freedom, and He couples 
the two. Consider first Ha- The most general quadratic form 
for the adatom Hamiltonian can be written 

1 

Ha = u'^j^„^,dl^s°'dm', (A9) 

a—Q,x,y,z m,m' — — 1 

where is the identity, s^'^'^ are Pauli matrices, and 
""mm' — ['"m' ml* '° Hcrmicity. Using the transfor- 
mation rules obtained above, one can show that symmetry 
strongly constrains the allowed w"^ terms. For example, 
in the a — sector, i?^/3 implies that uj^^ must van- 
ish except when m — m', while R^ further requires that 
""m m — "-m -m' (Timc-rcversal constrains these couplings 
no further) The a = x,y,z terms can be similarly con- 
strained, leading to the expression 

m=0,±l 

+ \/2A^„(4s^d-i + dls+di + h.c.) (AlO) 

provided in the main text. 

Let us elaborate on the physics here in greater detail. If 
the adatom were in isolation, continuous rotation symmetry 
would further restrict e|„j| to be independent of m and require 
that AgQ = Aso. In this case the Ago, A^^ terms above simply 
encode an isotropic L • S spin-orbit interaction that split the or- 
bitals into a lower j = 1/2 doublet and a j = 3/2 quadruplet 
that is higher in energy by 3 Ago. Crystal field effects aris- 
ing from the graphene environment lead to m-dependence of 
e\,n\ and allow A^^ to differ from Ago. Note also that we have 
assumed that the adatom's s-orbitals (and all other inner lev- 
els) are far in energy from the Dirac points and can be safely 
neglected. This is indeed justified by our supercell density 
functional theory (DFT) calculations for both indium and thal- 
lium, where the 5 s and 6s electrons contribute significantly to 
the density of states only at several eV below the Dirac points. 
We have further ignored four-fermion interaction terms in Ha- 



This, too, is justified for indium and thallium, since both ele- 
ments donate their outer p electron to graphene and form non- 
magnetic configurations. In a perturbative tight-binding pic- 
ture, electrons hopping onto the adatom from graphene there- 
fore lead to a singly occupied p orbital, for which Coulomb 
repulsion is expected to be unimportant. Finally, we note that 
Ha is even under R^ symmetry even though we did not en- 
force this. In other words, the fact that the adatom is displaced 
away from the graphene sheet leads to no additional terms in 
the Hamiltonian Ha compared to the (fictitious) case for an 
H-position adatom residing directly in the {x,y) plane. This 
fact is crucial for the irrelevance of induced Rashba terms for 
graphene, as we discuss further below. 

For the graphene-only part of the Hamiltonian, we posit that 
the usual nearest-neighbor hopping model Ht is modified pri- 
marily through a change in chemical potential for the six sites 
adjacent to the adatom: 

6 

Hg = Ht-SfiY,4,Cr,- (All) 

i=i 

Clearly Hg is invariant under not only Rj^/s, Rx, and T, but 
also Rz despite this symmetry being broken by the adatom. 
(Terms which violate R^ necessarily involve spin-orbit cou- 
pling, which we have justifiably neglected because it is ex- 
ceedingly weak.) 

Since the adatom generically modifies the bond lengths be- 
tween carbon atoms in its vicinity, the hopping amplitudes 
will be modulated as well near the adatom. According to our 
DFT calculations, however, perturbations to graphene's lattice 
structure are quite small (the carbon bond lengths in the imme- 
diate vicinity of either indium or thallium change by less than 
1 %). We thus expect this effect to be minor and have neglected 
it for simplicity. We have also neglected in Hg any spin-orbit 
terms which directly couple the carbon atoms, which is orders 
of magnitude weaker than the spin-orbit interaction indirectly 
mediated by the adatom. 

In contrast, the induced chemical potential J/i is not a weak 
effect — we find that our tight-binding simulations best repro- 
duce the DFT band structure when is on the order of 
leV. Physically, (5/i appears because the indium or thallium 
adatoms donate their outer p-orbital electron to the graphene 
sheet, leaving behind a net positive charge. Electrons in the 
graphene sheet thus prefer to conglomerate in its vicinity in 
order to screen the positively charged adatom. We note that 
of course there is no fundamental reason why the induced 
chemical potential should be confined to only the nearest six 
sites to the adatom. Our DFT simulations demonstrate, how- 
ever, that the induced charge modulation occurs very locally 
around the adatom, so that this is expected to be a good ap- 
proximation. As an illustration, consider the geometry with 
one thallium adatom in a 4 x 4 supercell. Figure |5] displays 
the induced charge density Ap = /9(graphcnc + thallium) — 
/9(graphcnG) — ^(thallium) relative to the charge densities 
obtained when graphene and thallium decouple completely. 
Here the yellow/blue regions correspond to isosurfaces with 
Ap — ±0.01eV/A'^. Clearly the charge modulation occurs 
predominantly within the first two 'rings' of carbon sites sur- 



9 



(a) 




(b) 




FIG. 5: Adatom-induced charge modulation, (a) Top and (b) side 
views of the charge density induced by a thallium adatom in a 4 x 4 
supercell. Yellow/blue surfaces correspond to a positive/negative in- 
duced charge density (relative to the case where graphene and thal- 
lium decouple completely). The charge modulations relax rather 
rapidly away from the adatom and occur mostly within the two in- 
nermost 'rings' of carbon sites near the adatom. 



rounding the adatom. 

Finally, let us discuss the Hamiltonian He that allows elec- 
trons to tunnel between the adatom and graphene sheet. For 
simplicity, we only allow for tunneling events which couple 
the adatom to its six nearest carbon atoms, and assume (be- 
cause of carbon's weak intrinsic spin-orbit coupling) that such 
processes are spin-independent. As in the main text it is useful 
to introduce linear combinations 
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1 



^(j-i). 



(A12) 



that, like rf„ 



under R, 



/3- 



, carry angular momentum quantum number m 
It follows from the symmetry transformation 



rules for obtained above that Cma transforms as follows: 



^7r/3 '■ 


Cm 




(A13) 


Rx ■ 


Cm 




(A14) 


T: 


Cm 




(A15) 


Rz ■■ 


Cm 




(A16) 



Using the first three symmetries along with the corresponding 
transformation rules for dm, it is straightforward to show that 



the spin-independent hoppings of interest must take the form 

Hc = - {t\m\Cldm + h.c.), (A17) 

m=0,±l 

with to real and ti pure imaginary. Under R^ symmetry, the 
to hopping is even while the ti hopping terms are odd. This 
is easy to understand physically by considering the fictitious 
case where the adatom resides exactly in the graphene plane 
so that Rz symmetry is present. In this case the ti terms reflect 
hopping between the adatom's p^.y orbitals and the carbon pz 
orbitals, which are orthogonal in this pathological limit. That 
is, tunneling between graphene's pz orbitals and the adatom's 
TO = ±1 orbitals arises only when Rz is lifted. In contrast, the 
adatom's pz orbitals overlap nontrivially with carbon's even 
when Rz is preserved. 

Remarkably, the only terms in H that do not respect Rz 
are precisely these ti couplings. To see what kind of spin- 
orbit processes the adatom can mediate, let us now consider 
the coupling between graphene and the adatom from a pertur- 
bative perspective. An electron from the graphene sheet can 
tunnel onto the adatom via either io or i i . Due to the adatom's 
strong spin-orbit coupling, the electron need not return to the 
graphene sheet via the same tunneling process by which it ar- 
rived. For example, an electron can tunnel via the hopping, 
flip its spin, and then return to the graphene sheet via ti. Con- 
sequently, at second order in perturbation theory, the adatom 
will mediate effective couplings between the carbon atoms 
that are proportional to tg, t\, and ^o^i- The electron's S^' spin 
component is conserved in the former two cases, but not the 
third. Furthermore, since the to coupling in He is even under 
Rz while the ti coupling is odd, it is only the induced terms 
for graphene proportional to t^ti that violate Rz- It follows 
that such spin-flip processes mediate Rashba spin-orbit cou- 
pling, while the t\ events mediate intrinsic spin-orbit interac- 
tions. (The tp processes do not generate spin-orbit coupling.) 
It is worth noting as an aside that our ability to classify such 
processes in this simple way relies on the fact that all terms in 
Ha are even under Rz- If Ha also involved terms that violated 
Rz, then any second-order tunneling event would generically 
be capable of mediating Rashba spin-orbit coupling. 

As discussed in the main text, to generate a quantum spin 
Hall (QSH) state in graphene, the physics near the Dirac 
points must be dominated by intrinsic rather than Rashba spin- 
orbit interactions. An important feature of our proposal is that 
Rashba coupling mediated by an adatom with a partially filled 
p shell is always subdominant to the induced intrinsic spin- 
orbit terms. This follows because electrons with momentum 
near the Dirac points hop very ineffectively onto the adatom's 
TO = orbital. Previously we discussed how this effect 
arises because electrons destructively interfere when tunnel- 
ing onto this orbital. Here we provide a complementary per- 
spective using the language of low-energy continuum Dirac 
fields. Focusing on physics near the Dirac points, one can ex- 
press the lattice fermion operators for the A and B sublattices 
of graphene as follows. 



(A18) 
(A19) 
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Here tpR/La slowly varying, two-component fields that 
describe low-energy excitations in the vicinity of the Dirac 
points at momenta ±Q = ± ^^^ x (a is the separation be- 
tween neighboring carbon atoms). Let us use this decomposi- 
tion to rewrite the operator Cm=o that appears in the to tun- 
neling. Neglecting pieces involving derivatives of tjjR/La^ we 
obtain 



1 

71 



-iQ ri 



[e "-■^ +e "--^ +e "-''^)V'l,b 



-iQrs 



(A20) 



Remarkably, all terms in parenthesis above vanish identically. 
The leading terms in the low-energy expansion of C„i=q thus 
involve derivatives and become, consequently, increasingly 
unimportant as one approaches the Dirac points. 

This implies that the Rashba coupling mediated by the 
adatom — which necessarily involves tunneling via to — is 
highly ineffective at influencing the physics at the Dirac 
points. In contrast, C„i=±i admits a non-trivial expansion 
in terms of iI^r/lo 6ven when terms involving derivatives 
are dropped. It follows that electrons with momenta near 
the Dirac points can effectively tunnel onto and off of the 
adatom via ti, and the intrinsic spin-orbit coupling mediated 
by these hoppings thus always dominates over the induced 
Rashba terms. As noted in the main text, this is not neces- 
sarily the case for adatoms with partially filled d or f shells, 
because there are then additional orbitals onto which electrons 
from graphene can tunnel. For instance, here electrons can 
tunnel onto the adatom's m — ±2 orbitals via a tunneling 
t2', such processes, unlike to tunneling, are not suppressed at 
the Dirac points. Effective couplings between carbon atoms 
proportional to tit2 are odd under R^, involve spin flips, and 
correspond to Rashba spin-orbit interactions that may signifi- 
cantly affect the low-energy physics at the Dirac points. The 
importance of such terms compared to the induced intrinsic 
spin-orbit interactions then depends on details of the Hamil- 
tonian. If, say, the additional orbitals hybridize very weakly 
with graphene or reside far in energy from the Dirac points 
due to crystal field effects, then d- or /-shell adatoms may 
still be suitable for stabilizing a QSH state in graphene. 



Appendix B: Density functional theory 

All DFT calculations were carri ed ou t with the Vienna 
ab-initio simulation package (VASPj^^EI at the level of the 
generalized-gradient approximation (GGaJ^. We used the 
projector augmented wave (PAW) method for the description 
of the core-valence interactiorF§E2l a vacuum space of 15 A 
was employed in all calculations, along with a 15 x 15 fc-grid 
mesh^" for integrals in the two-dimensional Brillouin zone. 
The energy cutoff of the plane wave expansion was set to 
500eV. 



For pure graphene, the optimized lattice constant is 2.468 A, 
consistent with the experimental value of 2.46A. When ad- 
dressing the preferred location (H, B, or T) of different ele- 
ments on graphene, a 4 x 4 supercell was employed with one 
adatom per cell. Recall that we have two minimal criteria for 
stabilizing a QSH state in graphene: the adatoms should (i) 
prefer to occupy the H position and (ii) not form a net spin 
moment. Among the elements we searched in this work. In, 
I, La, Hf, Au, Hg and Tl satisfy the first criterion, while the 
remaining elements prefer either the T or B site. However, 
I, La, Hf and Au form magnetic moments and thus fail the 
second criterion, while the s and d bands of Hg reside far be- 
low graphene's Fermi level, implying a negligible interaction 
between the two. The only two remaining elements — indium 
and thallium — satisfy both criteria and enter the H site non- 
magnetically. 

Including spin-orbit coupling, the binding energies [Et, = 
i?(adatom + graphene) — i?(graphene) — £'(atom)] of in- 
dium and thallium on graphene at the H site are -0.525 and 
-0.133 eV, respectively. The total energy per site when in- 
dium resides at the T and B sites is higher by 87 and 80meV, 
respectively, than when indium resides at the H site. This is 
qualitatively consistent with the results of Ref. 20 where spin- 
orbit coupling was neglected. The conclusion is the same for 
thallium, though the preference for the H position is weaker; 
here the T and B sites are higher in energy by 32 and 22meV. 
At low temperatures compared to these energy differences in- 
dium and thallium are expected to occupy only the H position, 
with indium-carbon and thallium-carbon bond lengths of 2.83 
and 2.90A as determined by GGA calculations. We also cal- 
culated atomic structures of indium and thallium on graphene 
within the local density approximation (LDA). The indium- 
carbon and thallium-carbon bonds are shortened by only 0.09 
and 0.06A, respectively, compared with the GGA results. It 
is thus reasonable to conclude that both GGA and LDA can 
describe the structural and electronic properties well. 



Appendix C: Adiabatic continuity of the adatom and 
Kane-Mele models 



We now turn to a periodic system with one indium or 
thallium adatom in a 4 x 4 supercell. This problem was 
explored extensively in the main text using both DFT and 
tight-binding simulations. Our first-principles calculations 
demonstrated that, with either element, the Dirac points re- 
main massless in the non-relativistic limit but acquire a sub- 
stantial gap upon including spin-orbit coupling. We further 
showed that the band structure and local density of states at the 
adatom site could be well-accounted for by an effective tight- 
binding model HixA- This model included first-, second-, and 
third-neighbor hopping amongst carbon atoms in the graphene 
sheet, as well as the additional couplings discussed in the pre- 
vious section around each adatom. Note that inclusion of 
these further-neighbor tunnelings is necessary for obtaining 
reasonable quantitative agreement with first-principles band 
structure at energies away from the Dirac points, even for 
pure graphene^'. Since the adatoms' local density of states 
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FIG. 6: Adiabatic continuity. Energy gap at the Dirac point as the 
Hamiltonian adiabatically deforms from the periodic adatom Hamil- 
tonian -ff4x4 at A = to the Kane-Mele model at A = 1. The solid 
blue and dashed red lines correspond to H4X4 evaluated with pa- 
rameters appropriate for indium and thallium, respectively. In both 
cases the gap remains finite as A varies from to 1, indicating that 
the Kane-Mele model and adatom Hamiltonian are adiabatically con- 
nected and thus support the same quantum spin Hall phase. 



is concentrated at energies greater than leV above the Dirac 
points (see Fig. 2 of the main text), one needs these addi- 
tional hoppings to accurately model their hybridization with 
the graphene bands. (We stress that this is important for a 
quantitative description only. None of the qualitative features 
discussed in this paper depend on these fine details. For in- 
stance, with only nearest-neighbor hopping a gap still opens 
at the Dirac points, though the local density of states for the 
adatoms' states is difficult to accurately fit in this case.) 

One advantage of this effective tight-binding model is that 
it enables one to readily verify that the insulating regime gen- 
erated by the spin-orbit-induced gap indeed corresponds to a 
QSH state. This can be perhaps most easily proven by demon- 
strating that iJ4x4 and the Kane-Mele model Hkm [see Eqs. 
(1) and (2) of the main text] can be smoothly deformed into 
one another without closing the Dirac point gap. To this end, 
consider the Hamiltonian H{\) = (1 — A)_ff4x4 + XHkm 
introduced in the main text that interpolates between these 
models as A varies from to 1. For H^m (which includes 
only nearest-neighbor hopping and uniform intrinsic spin- 
orbit coupling for graphene), we choose t = 2.7eV and 
Ago = — 0.0015eV. As Fig. |6] illustrates, the Dirac point gap 
computed for _ff (A) indeed remains finite for all A between 
to 1 with either indium (solid blue line) or thallium (dashed 
red line) parameters input into /f4x4- The QSH state known 
to be supported by Hkm and the insulating state stabilized by 
either type of adatom are, consequently, the same topological 
phase of matter. 



Appendix D: Graphene-only effective Hamiltonian 

In this section we outline the derivation of the graphene- 
only Hamiltonian, Eqs. (8-9) in the main text, from the tight 
binding model in Eqs. (4-7) describing the graphene sheet 
with H-position adatoms. We assume, for the purposes of this 



section only, that one adatom resides at each hexagonal pla- 
quette but disallow any direct hopping of electrons between 
adatom p orbitals. Under these assumptions we may exploit 
the translation symmetry of the problem and write the Hamil- 
tonian (4) in momentum space as H — J2k V'k^kV'k with 

V'k = (Ck,c^;c?ik,(^ok,rf-ik)^, where c^^^ remove elec- 
trons from the A/B graphene sublattices while c?mk remove 
electrons from the adatom orbitals. As before, spin is treated 
implicitly. The Hamiltonian matrix in this basis reads 



where 



Tk 

Tt 



hg T 



3=0 



(Dl) 



(D2) 



is the standard nearest-neighbor tight-binding Hamiltonian for 
graphene with 6j = (sin ^ , — cos unit vectors pointing 
from a site on sublattice A to its three neighbors on the B 
sublattice. Also, 



\ V2K,s^ 



describes the adatom p orbitals. Finally 




(D3) 



T = 



2 



(D4) 

represents transitions between graphene and the adatom or- 
bitals as described by He in Eq. (7). 

We now wish to 'integrate out' the adatom degrees of 
freedom and thus determine their effect on the electrons in 
graphene. The simplest way to accomplish this goal is to first 
perform a unitary transformation "H — > "H = e^^He^ with 
S chosen such that hg and ha in the transformed Hamiltonian 
are decoupled, i.e. T = 0. In this basis the adatom degrees 
of freedom can be integrated out trivially. Following the steps 
outlined in Ref.(9jwe find, to second order in T, 

hg = hg~^ [Th-^T^ + hgTh-^T^ + h.c] + O (T^) . 

(D5) 

This result is equivalent to treating Hamiltonian (Dl 1 pertur- 
batively to second order in T. 

We are interested in the effect of adatoms on the low- 
energy fermonic modes in graphene occuring in the vicinity 
of the Dirac points ±Q. In this limit, clearly, the second 
term in the brackets becomes negligible (since hg vanishes 
as k — > ±Q), so we therefore focus on the first term. Al- 
though it is possible to find an exact inverse of ha the result is 
cumbersome and tends to obscure the simple physics underly- 
ing the formation of the spin-orbit-induced gap in the system. 
To avoid this complication we make an additional assump- 
tion that spin-orbit splitting of the adatom orbitals is small 
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compared to the crystal field effects, |Aso|, jAg^j <C \em\, and 
write ha — ho + hi with ho — diag(ei, eo, ei). To first order 
in Ago and A^^ we then have h'^^ ~ h^^ ~ h^^hih^^ 



and 



h,-Th-'T^ 



Th-^hih-^T^. 



(D6) 



The adatom-induced correction to hg is seen to have a simple 
structure. The first term is spin-independent and mediates on- 
site potential terms and hoppings between sites in the hexagon 
surrounding an adatom. The second, spin-dependent term in 
turn produces the intrinsic and Rashba spin-orbit couplings. 



Various terms resulting from Eq. ( D6 1 can be explicitly 
worked out in a somewhat tedious but straightforward calcu- 
lation. The spin-orbit terms take the form 



soc _ / -^1 



Ri + D i?2 



rI Ri - D 



(D7) 



with 



\ti\' ' 



£1 = Ag„-^:^s^^sink- (Jj - Jj+i) (D8) 



and 



Ri = 



AL^?^ Y.[s X [dl - d,)l sink . {S, dl), 



3eoei 



R2 = -lAL 



/2to\ti\ 



3eoei 



If we now recall that Sq + Si + 82 — and identify 



Asn — A„ 



1^1 



Ar — A' 



V2to\ti\ 
3eoei 



(D9) 



we see that D and i?i 2 correspond respectively to the intrin- 
sic and hexagon Rashba spin-orbit couplings in Eq. (9). One 
can readily verify that R12 vanish at k = ±Q while D ap- 
proaches a nonzero value of 3\/3Aso, as expected on the basis 
of arguments presented above. 

Using the tight-binding parameters obtained in Sec. Ill of 
the main text we may estimate Ago — 23meV and A/j — 
58meV for thallium (for indium both are about factor of 3 
smaller). At 6% coverage, assuming uniform averaging of 
Ago over all plaquettes, this would imply a spin-orbit-induced 
gap of 0.06 X 6\/3Ago ~ 14meV, somewhat smaller that the 
21meV gap predicted by DFT at the same coverage. We at- 
tribute this discrepancy to the neglect of higher-order terms 
in Eq. (D61; since ^o.i are of similar magnitude as eo.i this 



expansion cannot be expected to yield quantitatively accu- 
rate values of the coupling constants. Nevertheless the pro- 
cedure outlined above is useful in that it illustrates how var- 
ious symmetry-allowed spin-orbit terms emerge in the effec- 
tive graphene-only model. It also confirms that although the 
Rashba coupling exceeds the intrinsic coupling by more than 
a factor of two, its special form, dictated by symmetry, forces 
it to vanish at the Dirac point and renders A^ irrelevant for the 
low-energy physics. 



Using the same methods one can also estimate various spin- 
independent terms induced by the adatoms. These include 
the on-site potential and additional hopping terms mentioned 
above. Since the latter are already present in pristine graphene 
and since adatom corrections are generally small compared to 
these bare values, we expect their effect on the low-energy 
physics to be minimal and do not discuss them further. 



Appendix E: Hexagonal Rashba and Disorder effects 

As mentioned in Sec. IV of the main text, the stability of the 
QSH phase was also tested against the addition of hexagonal 
Rashba coupling and residual disorder in the graphene sheet. 
We now discuss these results in greater detail. Throughout this 
section we will consider random adatom arrangements at cov- 
erage fii — 0.2, modeled using the graphene-only Hamilto- 
nian from Eq. (8) of the main text with (5/i = for simplicity. 
We first consider the influence of the hexagon Rashba term de- 
fined in Eq. (9). While it is true that its effects identically van- 
ish at the Dirac points, it nevertheless causes some splitting 
of the bands in their vicinity. As a result these Rashba terms 
can reduce somewhat the size of the mobility gap observed 
with either periodic or random adatoms, though to a far lesser 
extent than the traditional Rashba coupling. The splitting will 
depend on the ratio of Ago to A/j, and, more importantly, the 
absolute strength of Xji compared to t. 

In Figure [TJa) we show the result for a set of parameters 
chosen to illustrate this effect while avoiding the complica- 
tions due to the finite size of the system. Here Ago = 0.04< and 
Ai^ = cost, their ratio chosen to be similar to that estimated 



in Eq. ( D9 1 above. The mobility gap has clearly been some- 
what reduced, but a robust plateau remains. We remark that 
the standard Rashba term (with only nearest neighbor cou- 
pling) at the same coupling strength would completely destroy 
the plateau. For smaller values of Xr and Ago (but keeping 
their ratio constant) the effect of hexagonal Rashba coupling 
becomes even less pronounced although the finite size effects 
prevent us from obtaining reliable results at low adatom cov- 
erage in this regime. It is also worth emphasizing that one 
should not view Xji as reducing the gap compared to those 
found in our DFT and tight-binding results displayed in Fig. 2 
of the main text. Those simulations incorporated both the in- 
trinsic and hexagon Rashba spin-orbit couplings (among other 
couplings neglected for simplicity here), and the results al- 
ready reflect the influence of the latter. 

To study residual disorder unrelated to the adatoms (e.g., 
arising from impurities in the substrate), two models were 
used, both of which can be captured by an on-site potential 
term of the form specified in Eq. (8) of the main text. In the 
first model, we incorporate an uncorrected random chemical 
potential in the range [—A, A] on every lattice site. In the sec- 
ond, we ernploy a longer-range, correlated disorder potential 
of the forirP2l 
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FIG. 7: Hexagonal Rashba effect and substrate-induced disorder. Conductance G as a function of the Fermi energy Ef for a system 
including (a) hexagonal Rashba coupling, (b) uncorrelated disorder randomly distributed over the range [—A, A] and (c) correlated disorder 
as defined in Eq. l |El[ ) with p — 0.07 and d = 5. In all cases the system is of size W — SO and L = 40 with adatom coverage rii = 0.2 and 
nearest-neighbor hopping strength t = 2.7eV. For panel (a) the intrinsic spin-orbit strength is Aao = 0.04i; for panels (b) and (c), Aso = Q-U. 



where -/Vimp is the number of impurity sites, Vj is also uni- 
formly random over the range [—A, A], and d is the radius of 
the impurity potential. As in Ref. 52^, we characterize iVimp 
by the ratio p = Ni^^p/N, with N being the total number of 
carbon atoms in the sample. 

While we have not averaged over multiple disorder realiza- 
tions as in Fig. 3(b) of the main text, with only a single disor- 
der realization one can still see the general trend towards insta- 
bility of the QSH phase as A increases. Our simulations show 
that this state, despite already being stabilized by a disordered 
arrangement of adatoms, in fact exhibits remarkable resilience 
against additional uncoiTelated substrate-induced disorder in 
the graphene sheet. As Fig.|7jb) illustrates, in this case a very 



high value of A « i can be applied before the first signs of 
disorder affecting the topological phase appear, and even here 
a considerable conductance plateau survives. The topologi- 
cal phase is, however, more sensitive to long-range correlated 
disorder, which is likely more relevant for experiment when 
charged impurities from the substrate present the main disor- 
dering mechanism. Figure [TJc) illustrates that here the deteri- 
oration of the plateau is more rapid. At A = Q.2t w 0.54eV, 
with p = 0.07 and d = 5, there is still a remnant of the gapless 
edge states, but for stronger disorder the topological phase is 
destroyed. Notice that this value of A is comparable to the 
width of the plateau in the 'clean' case where the only disor- 
der source is the randomly distributed adatoms. 
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